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, We investigate the inversion of perturbation series and its resummation, and prove that it is related 

^ , to a recently developed parametric perturbation theory. Results for some illustrative examples show 

, that in some cases series reversion may improve the accuracy of the results. 

(N 
Ph. 

r-] ! I. INTRODUCTION 

» 

Perturbation theory yields the solution of a problem in the form of a power series of a properly chosen model 
(or dummy) parameter. When the convergence radius of this series is too small for the physical application, or 
^ ■ the series converges too slowly, it is customary to resort to a resummation method that produces an approximant 
that improves the results. There are many approaches for that purpose; among them we mention Pade approximants, 
On ■ Borel summation, algebraic approximants, continued fractions, nonlinear transformations of the expansion parameter, 

J: etc. Q, 0,3 . 



PACS numbers: 03.65.Ge 



^-H " The purpose of this paper is to investigate the application of a well-known mathematical method, named series 
, reversion or series inversion [J, to perturbation theory. In Sec. |lT] we outline the main ideas of the approach, 

^ in Sec. IIIII we show the connection between the method of series reversion and a recently developed parametric 

■ perturbation theory 

Has 

. In Sec. IIVI we compare the accuracy of resummation of the direct and inverse series for 

■ some illustrative examples. Finally, in Sec. |V|we draw some conclusions about the usefulness of the approach. 

II. INVERSION OF SERIES 

Suppose that we are trying to estimate accurate values of an unknown function E{g) from a few available coefficients 
of its power series 

oo 
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When the convergence radius of this series is too small for the physical application, or the series converges too slowly, 
it is custommary to resort to a resummation method that produces an approximant E w A^^^ (g) from the partial 
sum of order N: Sn = Eq + Eig + . . . + E^g 

It is always possible to invert the series ([Ij and obtain g in terms of AE = E — Eq: 

^ AE E-iAE^ (2E^ - E1E3) AE^ 

9 = AeY^G,AE^^ ^^^ ^ ; (2) 

If we apply a resummation method to this series we obtain an approximation of the form g w (^^)- This strategy 
may give more accurate results than the former one if the radius of convergence or the region of utility of the series 
([2|) is greater than that of the direct expansion ([T]). However, even in this favorable situation we are paying the price 
of having the inverse of the desired function. 

In order to get some motivation for the use of the inverse series, consider the function 

E = yiT^ = 1 + _ ig2 ^ . . . (3) 

that is real for all g > —1. The Taylor series about 5 = converges only for j^j < 1 because of the branch point 
at g = —1. However, the inverse series g = 2AE + AE^ converges for all AE suggesting that in some cases it is 
convenient to use the latter instead of the former. 



III. PARAMETRIC PERTURBATION THEORY 

In what follows we show that the inversion of series is related to a recently developed parametric perturbation 
theory 0,3,0]. If we define the parameter p = AE/Ei then 

g^p^GjEi p'=p-—- + -2 + ... (4) 

j=o ^ 1 

Application of a resummation method gives an expression for g in terms of the parameter p: 

g^A[^\E,p) (5) 

which toghether with 

E = Eo+ Eip (6) 



define an approximate parametric representation of E(g)- These equations are the basis of the so called parametric 
perturbation theory developed recently by Amore 0, from a principle of absolute simplicity. 

More precisely, in parametric perturbation theory one introduces an approximant g ~ A^^^ (p) into the direct 
expansion ([1]) and choose the approximant coefficients so that E = bo + bip + Since the approximant is 

constructed so that A^'^^p) w p for p 1 0, 0, 0] then we realize that Bq = Eq and bi — Ei, and conclude that 
parametric perturbation theory is basically the inversion of the perturbation series. 

In order to see the connection clearer consider the trivial example ^ again. If we substitute p + C2p^ + . . . + cmP^ 
for g into the series ([3]) and choose the coefficients Cj in order to remove the terms of order greater that one, we obtain 
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C2 = 1/4 and Cj = 0, j > 2. The resulting expressions E = 1 + p/2 and g = p + exactly agree with the inverse 
series shown above. 

In some cases one easily estimates a range of utility of the parametric representation of E{g). For example, suppose 
that we know that dE/dg < (> 0), then i?i < (> 0). In both cases the range of utility of the approximant 
g w A^^\p) is determined by the conditions p > and dg/dp > 0. 



IV. EXAMPLES 



As an illustrative example we consider the integral 

/•OO 

E{g) = / e-^'-9^*dx. (7) 
Jo 

The formal expansion 

, A _ 3g 105.g^ _ M65g^ 675675g4 _ 436486053^^ \ 

2 V 4 32 128 2048 8192 "') ^ ' 

is known to diverge for all values of g 2] which is reflected by the form of the expansion coefficients En = (— l)"r(2n + 
l/2)/[2(n!)]. We can rewrite the inverse expansion 

8AE 280A£;2 

^(^^)--37^ + ^^-- 

in terms of the parameter p = —S^/ttAE/S as 

35p2 35p3 i7675p4 263095^^ 
^(^^=^+^ + ^ + ^56 25^ + --- ^''^ 

In this case we have 

Eip) - f - ^ (11) 

The pair of equations for g(p) and E{p) are a parametric representation of the function E{g). 

For simplicity we restrict to series of order five with the purpose of illustrating some ways of constructing the 
approximants. In order to determine the range of applicability of the p-series we take into consideration that dg/dp > 
for all p > 0. Therefore, we assume that it may be resonable to use the parametric representation for all < p < pm, 
where pm is the smallest positive root of dg/dp = 0. In the case of the partial sum of order five pop we have 
Pm = 0.1659 that leads to = g{pm) — 0.2194. Table [J shows that the parametric representation embodied in 
equations (fTO|l and (fTT|) yields better results than the direct g-power series ([8|), at least for all g < gm- 

We can improve the accuracy of the results by means of resummation methods. For simplicity we try Fade 
approximants built from the ^-series and p-series of order five shown above. Since E{g oo) = Q we choose 

the [2/3] Fade approximant 

^ _ 20F(11546O139.g2 + 28532448.g + 1163200) 

^ 141105195.g3 + 5347881005^ + 117619392g + 4652800 ^ ' 

The [2/3] Fade approximant for the inverse series 



16p(219707p+ 15640) 
^ ~ 26152805p3 - 11137140p2 + 2420512p + 250240 ^ ' 
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presents a maximum at pm = 0.360 where = 0.607. Table [Til shows that the Pade approximant for the direct series 
is more accurate than the one for the inverse series. 
The [3/2] Pade approximant 

_ p (15904 + 318204p + 747223p2) 
^ ~ 15904 + 248624p - 375297p2 ^ ^ 

exhibits a pole at po = 0.721 and we can therefore use this parametric representation for all < 5 < oo when 
< p < po- The main limitation of this approach is that E{g — > 00) = Eq + -Eipo = 0.407. Table Hill shows that the 
inverse series is slightly more accurate for small g but leads to a wrong limit. On the other hand, the direct series 
approaches the correct limit too fast as 1/g. 

The approximants so far used do not take into consideration the asymptotic behavior E{g) ~ as g 00. In 

the case of the direct series we may calculate [N/N + l]{g) Pade approximants for E{g)'^ and then use [N/N + \\{gY/^ 
as a reasonable approximation; for example 

/ 116949g^ + 27216.g+ 1060 \ ^''^ 
2 ^^218277^3 + 19064752 + 303965+ 1060/ ' ^ ^ 

We do not know how to obtain reasonable approximants for the inverse series; for that reason we try Amore's approach 
g[p) = p[N/N + l]{p)^; for example, 

-1048576p(24078737127p2 + 4884321032p + 161747680)'^ 
^ ~ (240656732281p3 - 321691293064p2 - 75884668992p- 2587962880)^ ^ ' 

constructed from the p-series of order six. One realizes that this is the parametric perturbation approach proposed 

by Amore HQS 

for this particular case. In the neighbourhood of the pole p — pa closest to origin this approximant 
behaves as 5 « Koip — Po)^^ so that E{g) ^ Eq + Eip^ + EiK^/^ g^^^^ . We appreciate that this parametric 
representation tends to a wrong limit (because Eq + Eipo 7^ in general) and in a wrong way. Table [IV] shows that 
these approximans are more accurate than the previous ones, as expected, and that the direct series is clearly better 
than the inverse one. Notice that the parametric representation gives the wrong limit E{g — > 00) — 0.861. Besides, 
remember that we need perturbation coefficients of order 5 and 6 of the direct and reverse series, respectively, in the 
construction of such approximants. 

Following the same philosophy we have also tried [3/2] approximants for the reverse series 

_ p (583627586759p3 + 767002386296p2 + 44364247040p - 265333760)^ 
^~ 32768 (290961289397p2 + 5574551760p- 33166720)^ 

Table IVl shows that this parametric representation yields better results than the preceding one, but it also leads to a 
wrong limit: E{g 00) — 0.883. 
The polylogarithm function 

^^^(^) = T>)j/ 

appears in several fields of theoretical physics, for example, in the Bose-Einstein and Fermi-Dirac distributions [lol ]. 
The Taylor expansion of this function about z — yields the series 



Lz.(z) = ^- (19) 
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In this case the inverse series gives more accurate results than the direct one and, consequently, we draw the same 
conclusion regarding the Pade approximants constructed from them. Table I VII shows that the Pade approximant 
[5/6] (p) for the former series gives considerably more accurate results than the [6/7] (5) Pade approximant for the 
latter one. 

The formal Taylor series about g = for the function 



is divergent as shown by its coefficients £'„ = (— l)"n!. The asymptotic behavior at 5 3> 1 is given by 

+ + (21) 

9 9 

where 7 = 0.5772156649 ... is Euler's constant. In this case the inverse series is more accurate than the direct one, 
and exactly the same situation takes place for their corresponding Pade approximants. 

As another illustrative example consider the function E{g) defined hy g = Ee~^. The power series E{g) = 
9 + 9^ + 3.9^/2 + . . . converges for g < which is reflected by the form of the expansion coefficients £'„ = n"^^/n!, 
n > 1. On the other hand, the inverse series g — E — E"^ + . . . + (— l)"^^i?"/ (n — 1)! + . . . converges for all E and 
is therefore preferable for practical applications. However, it is not always true that the inverse series have greater 
convergence radius than the direct one; simply consider the function E{g) — ge^^ , where the role of the variables has 
been reversed. In this case parametric perturbation theory will perform poorer than standard perturbation theory. 



V. CONCLUSIONS 

In this paper we investigated the usefulness of the inversion of perturbation series and its resummation. Our simple 
examples show that in some cases it is convenient to resort to the inverse series, but in others the straightforward 
perturbation series is expected to provide better results. We have also proved that the recently proposed parametric 
perturbation theory , based on the principle of absolute simplicity, consists of a convenient modification 

and resummation of the inverse series. Therefore, in some cases this approach will not perform better than well- 
known resummation methods on the direct series. In addition to it, parametric perturbation theory exhibits two 
disadvantages: producing the inverse of the desired function, and the difficulty of deriving the correct asymptotic 
limit and behaviour when it is known. 



[1] C. M. Bender and S. A. Orszag, Advanced mathematical methods for scientists and engineers (McGraw-Hill, New York, 
1978). 

[2] G. A. Arteca, F. M. Fernandez, and E. A. Castro, Large order perturbation theory and summation methods in quantum 

mechanics (Springer, Berlin, Heidelberg, New York, London, Paris, Tokyo, Hong Kong, Barcelona, 1990). 
[3] F. M. Fernandez, Introduction to Perturbation Theory in Quantum Mechanics (CRC Press, Boca Raton, 2000). 
[4] L. V. Ahlfors, Complex Analysis (McGraw-Hill, New York, 1953). 

[5] P. M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York, 1953). 

[6] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Ninth ed. (Dover, New York, 1972). 



6 

[7] P. Amore, A new approach to resummation: Parametric Perturbation Theory. larXiv:0705. 11661 
[8] P. Amore, Going beyond perturbation theory: Parametric Perturbation Theory, a rXiv:0705. 11681 
[9] P. Amore, Phys. Rev. D 76, 076001 (2007). 
[10] D. A. McQuarrie, Statistical Mechanics (Harper & Row, Publishers, New York, 1976). 



TABLE I: Series expansions of order five for the integral (O 
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P 


-cja.i|o|i 


T7„ fli Ah 
Jiq.lliyjl 


Exact 


0.04 


0.227268254 


0.8630223905 


0.8632172917 


0.8632281022 


0.08 


0.2211979138 


0.8358839701 


0.8447690498 


0.8449941504 


0.12 


0.213865426 


0.7500155805 


0.8285524745 


0.8296883134 


0.16 


0.2043179764 


0.4525374722 


0.7504228564 


0.8164228001 


0.2 


0.1891655735 


-0.3655376234 


0.7604942069 


0.8046805576 



TABLE II: Pade approximants [2/3]for the integral © 

g p Eq.llll) Eq.lig Exact 

0.1 0.07419851329 0.836884189 0.8369093852 0.8370429277 

0.2 0.124755604 0.803286224 0.8033055939 0.8046805576 

0.3 0.166889096 0.7760554974 0.7753007175 0.7800434542 

0.4 0.2069906611 0.7523625938 0.7486464025 0.7600358198 

0.5 0.2515059286 0.7310124346 0.719058431 0.7431554088 

0.6 0.392980645 0.7113938594 0.6250244038 0.7285439429 



TABLE III: Pade approximant [2/3] for the (/-series and [3/2] for the inverse series of the integral 

g p Eq.jn} Eq.llll) Exact 

0.1 0.07404520814 0.836884189 0.8370112825 0.8370429277 

1 0.3079042644 0.6446831837 0.6815721382 0.6842134278 

10 0.6043806547 0.2142628036 0.4845131183 0.4609804743 

100 0.7061454902 0.02801405374 0.4168730654 0.2772884009 

1000 0.7196727352 0.0028904009730.4078819088 0.1594808649 

10000 0.7210732433 0.0002899613860.4069510328 0.09033502245 



TABLE IV: Improved Pade approximants [2/3] for the integral ([7| 

g p Eq.llS} Eq.lHl) Exact 

0.1 0.07356402668 0.8369445716 0.8373311095 0.8370429277 

1 0.2892446121 0.6715809835 0.6939746529 0.6842134278 

10 0.5828537448 0.4198251659 0.4988214137 0.4609804743 

100 0.03723484095 0.2393860646 0.8614780364 0.2772884009 

1000 0.03737319161 0.1348101326 0.8613860789 0.1594808649 

10000 0.03737319161 0.07582023369 0.8613860789 0.09033502245 
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TABLE V: Improved Pade approximant [3/2] for the inverse series of the integral ([7| 



Eq.dUl) 



0.1 0.07404345742 0.8370124462 

1 0.3170281954 0.6755077332 

10 0.8091513837 0.3484081181 

100 0.00545012921 0.882604387 

1000 0.00545012921 0.882604387 

10000 0.00545012921 0.882604387 



TABLE VL Values of Li3/2(ff) for g close to 1 obtained from Pade approximants on the direct and inverse series. Exact figures 
are underlined for comparison. 
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[6/7] (3) 


[5/6] (p) 


Exact 


0.999999 


2.380740506 


2.608744256 


2.608831900 


0.99999 


2.380591082 


2.601153011 


2.601179942 


0.9999 


2.379099267 


2.577063920 


2.577071427 


0.999 


2.364418183 


2.501706883 


2.501708465 


0.99 


2.237103024 


2.271659944 


2.271660077 


0.9 


1.614336255 


1.614438528 


1.614438529 



